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Abstract. It has been suggested that under certain approximations the Sivers effect can be described 
in terms of factorization of final state interactions and a spatial distortion of impact parameter space 
parton distribution; that is a convolution of the so-called lensing function and the impact parameter 
GPD E. In this approach the lensing function is calculated in a non-perturbative eikonal model. 
This enables a comparison between the a priori distinct Sivers function and the GPD E which goes 
beyond the discussion of overall signs. 

Keywords: Transverse Momentum, Gluonic Poles, TMDs, GPDs 
PACS: 11.15.Tk, 12.40.-y 

"T-odd" transverse momentum-dependent (TMD) parton distribution functions 
(PDFs) have gained considerable attention in recent years. Prominent examples, and 
the subject of extensive study, are the Sivers and Boer-Mulders functions 01 01 • In the 
factorized picture of semi-inclusive deep inelastic scattering Hi 0] at small transverse 
momenta Pr ~ k± « a/(2^ me Sivers effect describes a transverse target spin asym- 
metry (TSSA) through the "naive" T-odd structure, Af(x,k T ) ~ S T ■ (Px^)/^!,^). 
Dynamically, "T-odd"-PDFs emerge from the gauge link structure of the multi-parton 
quark and/or gluon correlation functions 10,0,01 which describe initial/final- state inter- 
actions (ISI/FSI) of the active parton via soft gluon exchanges with the target remnants. 
Though these interactions are non-perturbative recent phenomenological calculations 
approximate the FSIs by perturbative one-gluon exchange lf8L [9L ITol. ITll IT2L IT3I. [l4l IT5h . 



We improve this approximation [|16ll by applying non-perturbative eikonal meth- 



mp 

ods It 1 7L 1 1 8H to calculate higher-order soft gluon contributions. This approach provides 
insight into the picture of the Sivers effect proposed by Burkardt in the mixed phase 
space {x, bj} description of the nucleon [T^, 2(J. We consider the first kj moment of the 



Sivers function written as a gluonic pole matrix element in impact parameter space pll] 

(4«> = jd 2 b T j ^ r e" p+z "(P+,0 r ;5r|?(zi)y + [zi;z2]/ ! '(z 2 )9(z2)| J P + ,0r;5r) (1) 

where (k l T (x)} = Me('S J T f^ l \x). The operator / originates from the time -reversal be- 
havior of the FSIs implemented by the gauge link operator [a ; b] which denotes a straight 
Wilson line connecting the two locations a and b. It is a function of the gauge-link 

and the gluonic field strength tensor, l l {^ri) = / dy~ [^n 7 y~n]gF +l (y~n) \y~n, %-n] 

where = 1/^(1,0,0,-1) is the light-like vector, andz 1/2 = (=F^,^r, ±^75 We 




FIGURE 1. Left: TMD amplitude W(P,k;S) including FSIs for re-scattered eikonalized quark and the 
spectator. The final state interactions are described by a non-perturbative amputated scattering amplitude 
M. Right: FSI in the relativistic eikonal model where quark and spectator are eikonalized. Only ladder 
diagrams are taken into account. 



find that calculating the operator / in an eikonal approximation, including multiple gluon 
exchanges, the integrand of Eq. (OQ) factorizes into a distortion of the transverse space 
parton distribution and the FSIs resulting in the relation Q20Q 

-M 2 4s J T f^ 1 \x) ^ J d 2 b T J\xM)4vA dp T (#(x,b 2 T )). (2) 

Here J* represents the "chromodynamic lensing function". While such a factorization 
doesn't hold in general [22, 23] it is fulfilled in lowest order contributions of simple 
field-theoretical models [|24L |21|] . In particular, contributions beyond the lowest order, 
where the quark fields and the operator / "interact" (they are interacting Heisenberg 
operators) leads to a breakdown of the relation which we write as /' = e ^ , '(x,i>r)l+b.t., 
where b. t. denote the "breaking terms". We estimate the size of these breaking terms by 
deriving the lensing function in a relativistic eikonal model dial . 

If the FSIs are modeled by a one-gluon exchange, the lensing function can be unam- 
biguously identified in a diquark spectator model [24]. For this reason we choose to stay 
in the diquark picture but describe the FSIs by a non-perturbative (re-)scattering at the 
amplitude level as in Fig. (0Q). We perform a calculation of M and in turn the lensing 
function in a U(l) theory, and a SU (2)-non-Abelian theory in an attempt to study the ef- 
fects of color on the relation Eq. (O. Thus, we do not take into account gluon exchanges 
between the quark line and the soft blob since they give rise to interactions between the 
quark fields and the operator / in Eq. (OQ) which lead to a breakdown of the relation, 
Eq. ©. Self-interactions and real gluon emission between quark and diquark lines are 
accounted for through masses and vertex form factors in the spectator model. Applying 
the usual Feynman rules, and expressing the FSIs through the amputated quark-diquark 
scattering amplitude M leads to the expression 

AWYPttt / A iSN{{P-q) 2 )W-4 + m q )u{P^)]M{ q ,P-k) 

A ' ' J J (27t) 4 [n-(P-k-q) + i0}[(P-q) 2 -ml + i0][q 2 -m2 + i0} ' K) 

where i(n - (P — k — q) + /0) -1 represents the eikonal propagator, the double-line in Fig.[T] 
and AW — W — W° where W° denotes the contribution without final-state interactions. 
Performing the contour-integration on q~ we pickup poles only from the denominators 
in Eq. © where we assume that M does not contain poles in q . This is not generally true 
even in the one-gluon exchange approximation |25L 3h where poles exist for higher twist 



T-odd TMDs for a scalar spectator and also at leading twist for an axial vector spectator. 



However, these contributions give rise to light-cone divergences ll25ll which would lead 
to a breakdown of Eq. ©. Next considering the integration on q the eikonal propagator 
can be split into a real and imaginary part. While the imaginary part of the eikonal 
propagator forces the diquark to be on-shell and restores the relation, Eq. (O the real 
part leads to contributions beyond the relation and has to be attributed to the breaking 
terms. The resulting kinematics are consistent with the picture that quark and diquark 
move quasi-collinearly with respect to the nucleon and that the FSIs are dominated by 
"small" transverse momenta gluon exchanges. Now as input to Eq. © we calculate the 
GPD £(x,0, —A 2 ) for a scalar diquark spectator with a multi-pole form factor 



(27t) 3 J [(fc r -i(l-x)A r )2+A 2 ]"[(2r+5(l-x)A r )2+A 2 ]' 

where A, m q , m s and M are "cutoff", quark, spectator, and nucleon masses respectively 
and A 2 = xm 2 — x(l — x)M 2 + (1 — jc)A 2 . We transform to impact parameter space, 
S{x^bj) = J j^ I e^' Aj " bT E(x, 0, — A 2 ) and we express the lensing function J? in terms 
of the real and imaginary part of the scattering amplitude M, 

(5) 

where p and q are transverse momentum integration variables, and calculate M in the 
eikonal approximation (details can be found in lllql ) 

M eik (x, \q\) = Atz{1 -x)P + j~dzzJ (z\q\) (e^ - l) , (6) 

where % (z) is the eikonal phase in configuration space 

/OO /'CO 
da dp^(z+an-pn), (7) 
-00 J — 00 

n = l/v2( 1,0,0, 1), e q e s = —4na and 2$ P a(z) are the quark-gluon coupling and gluon 
propagator respectively. Eq. © is valid in an U(l) Abelian gauge theory. Using Eq. © 
for M eik and transforming to impact parameter space the lensing function 

<)'Wr) = ^-X'(U) (l +cos Z (§i)) . (8) 

An attempt to build in color into the eikonal formalism has been provided in ||18ll . Using 
these results the lensing function in an SU(2) gauge theory becomes 



(9) 



In principle, SU(3) color can be implemented numerically into the formalism, however 
it is non-trivial to derive an analytical expression similar to Eqs. d8]) and ©. The eikonal 



FIGURE 2. Left: Lensing function for the perturbative U(l) and non-Abelian SU(2) cases. Note 
that a value of a = 0.3 in the perturbative result is arbitrarily chosen. Center: Valence u-quark GPD 
E^ u '(x,0, —Aj) at x = 0.3 in spectator model iflql and from the parameterization [ 28] . Right: First moment 
of the u-quark Sivers function from Eq. (O for the two models Il28lll6ll of input GPD. 



phase in (Q contains two characteristic quantities of this formalism, the coupling a 
and the gauge propagator 3>~ x . They represent the exchanged soft gluons between 
quark and spectator. In order to numerically estimate the lensing function and in turn 
the Sivers function we need to implement the infrared behavior of the gluon and the 
running coupling in the non-perturbative regime where we infer that the soft gluon 
transverse momentum defines the scale at which the coupling is evaluated. These two 
quantities have been extensively studied in the infrared limit in the Dyson-Schwinger 
framework 112611 and in lattice QCD Il27n . We use calculations of these quantities from 
Dyson-Schwinger equations [26] where both a s and 3>^ x are defined in the infrared 
limit. With these inputs we plot the lensing function versus |Z?r| for the perturbative, 
U(l), and SU(2) cases in Fig. [2l We find that for all cases the lensing function is 
"attractive" as suggested in the general arguments of Burkardt |20(]. In order to study 
the dependence of the relation © on the GPD E, we compare the diquark model result, 
Eq. © to another phenomenological model investigated in Ii28ll . In Fig. [2] we plot the 
first moment of the Sivers function for the various cases which range from 0.04 — 0.08. 
This is in reasonable agreement with the extractions from the Torino and Bochum groups 
is PS 0.05 0, [30]. In this calculation we have shown that Eq. © is indeed preserved 
when higher order diagrams contribute to the lensing function under conditions of the 
eikonal and ladder approximation. We have also provided an improvement to the one- 
gluon exchange approximation in calculating T-odd TMDs in the spectator framework. 
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